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Propellant slosh is a potential source of disturbance critical to the stability of space vehicles. 
The slosh dynamics are typically represented by a mechanical model of a spring mass damper. 
This mechanical model is then included in the equation of motion of the entire vehicle for 
Guidance, Navigation and Control analysis. Our previous effort has demonstrated the 
soundness of a CFD approach in modeling the detailed fluid dynamics of tank slosh and the 
excellent accuracy in extracting mechanical properties (slosh natural frequency, slosh mass, 
and slosh mass center coordinates). For a practical partially-filled smooth wall propellant 
tank with a diameter of 1 meter, the damping ratio is as low as 0.0005 (or 0.05%). To 
accurately predict this very low damping value is a challenge for any CFD tool, as one must 
resolve a thin boundary layer near the wall and must minimize numerical damping. This 
work extends our previous effort to extract this challenging parameter from first principles: 
slosh damping for smooth wall and for ring baffle. First the experimental data correlated into 
the industry standard for smooth wall were used as the baseline validation. It is demonstrated 
that with proper grid resolution, CFD can indeed accurately predict low damping values from 
smooth walls for different tank sizes. The damping due to ring baffles at different depths from 
the free surface and for different sizes of tank was then simulated, and fairly good agreement 
with experimental correlation was observed. The study demonstrates that CFD technology 
can be applied to the design of future propellant tanks with complex configurations and with 
smooth walls or multiple baffles, where previous experimental data is not available. 


I. Introduction 

Propellant slosh is a potential source of disturbance critical to the stability of space vehicles. The slosh dynamics 
are typically represented by a mechanical model of a spring mass damper. This mechanical model is then included in 
the equation of motion of the entire vehicle for Guidance, Navigation and Control analysis. The typical parameters 
required by the mechanical model include natural frequency of the slosh, slosh mass, slosh mass center location, and 
the critical damping ratio. During the 1960’s US space program, these parameters were either computed from an 
analytical solution for a simple geometry or by experimental testing of sub-scale configurations. Our previous work 
[1] has demonstrated the soundness of a CFD approach in modeling the detailed fluid dynamics of tank slosh and 
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has shown excellent accuracy in extracting the mechanical properties for different tank configurations as a functiuon 
of fill level. The verification and validation studies included a straight cylinder against an analytical solution, and 
sub-scale Centaur LOX and LH2 tanks with and without baffles against experimental results for the slosh frequency, 
slosh mass and mass center. The study shows that CFD technology can provide accurate mechanical parameters for 
any tank configuration and is especially valuable to the future design of propellant tanks, as there is no previous 
experimental data available for the same size and configuration as the current flight designs. 

Since the liquid oscillatory frequency may nearly coincide with either the fundamental elastic body bending 
frequency or the dynamic control frequency of the vehicle at some time during the powered phase of the flight, the 
slosh forces could interact with the structure or control system. This could cause a failure of structural components 
within the vehicle or excessive deviation from its planned flight path [2], It is therefore necessary to consider means 
to provide adequate damping of the liquid motion and slosh forces and to develop methods for accounting for 
damping in the analyses of vehicle performance. Determination of slosh damping in a given tank configuration is a 
very challenging task. First, an analytical solution does not currently exist for the slosh damping due to high 
nonlinearity of the problem. While slosh frequency can be computed using linear potential theory, the damping 
physics involves the vorticity dissipation which requires full solution of the nonlinear Navier-Stokes equations. 
Previous investigations and knowledge of damping characteristics were all carried out by extensive experimental 
studies. Previously, four extensive experimental investigations have been carried out on viscous damping in a 
circular cylinder [3-6], and the damping values have been correlated to a functional form of: 


y= Cx/Re 


(l) 


where Re is a dimensionless parameter analogous to an inverse Reynolds number [7]: 


Re = 


/gR* 


( 2 ) 


and C is a constant, y is the damping ratio, or the critical damping ratio of the amplitude of the free surface 
oscillation, R is the tank radius, g is the gravitational acceleration, and v is the kinematic viscosity of the liquid. 
Mikishev and Dorozhkin [6] proposed the following correlation from their tests [7]: 
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Where h is the liquid depth. For large depth of h/R > 1 .0, the above equation may be approximated by: 

^=0.79x/Re 


(4) 


A similarly extensive but independent study by Stephens et al. [5] found a slightly different correlation: 
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When the liquid depth is large, equation (5) reduces to: 

r= 0.83VRe 


( 6 ) 
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The above correlations have become the industry standard methodology to compute slosh damping value. For 
example, NASA Engineering and Safety Center (NESC) performed a comparison of the Orion Service Module slosh 
damping values with the Ares-1 launch vehicle slosh damping values and legacy slosh damping values from the 
Space Shuttle Program and the Apollo Saturn launch vehicle [8]. The empirically computed 1-g minimum slosh 
mode damping value from the above correlation for the Orion nitrogen tetroxide (NTO) oxidizer tank was found to 
be 0.0003 or 0.03%. Likewise the empirically computed 1-g minimum slosh mode damping value for the Orion 
monomethyl hydrazine (MMH) fuel tank was found to be 0.0006 or 0.06%. The NESC assessment team verified 
these empirical 1-g calculations of damping values. 

It should be noticed that the above correlations are only for right cylinders and applicability to curved bottom 
tanks is questionable. With advancement of CFD technology, it is possible to predict slosh damping directly from 
first principles. As one may realize the critical damping ratio for practical partially-filled propellant tank is as low as 
0.03% to 0.06%, this gives rise to a great challenge for any numerical prediction using a CFD tool, as the numerical 
damping used to stabilize the solution could be higher than the viscous damping. Another challenge is the 
requirement to resolve thin boundary layers near the walls. 

The objective of this effort is to improve our understanding of the physics behind slosh damping and to validate 
CFD extracted damping against experimental data using the same parameters. The study will take a fundamentally 
sound approach first with validations against experiments for the smooth wall cylindrical tank. High-order numerical 
schemes will be applied using a technique developed to estimate and reduce/remove the numerical damping from 
the solution. With the validated CFD model, we will then study the damping in the presence of a flat ring baffle 
which is a commonly used as means of slosh suppression. 


II. Computational Modeling Tool 

The computational software used to study the tank vertical sloshing phenomenon is the commercially available 
CFD-ACE+ program, which was originally developed by CFD Research Corporation (CFDRC), and is currently 
owned and distributed by ESI [8], CFD-ACE+ is a multi-physics and multi-disciplinary simulation tool, and is 
especially suited for liquid slosh modeling. 

CFD-ACE+ solves the Navier-Stokes equations in a Lagrangian-Eulerian frame. The continuity and momentum 
equation can be generally written as: 

~ j pdV + { (p(v- v g )-ds) = 0 (7) 

dt V s 

— | p4>d V + J p(f>(v — v g ) - ds = J q • ds + J S^dV (8) 

dt V s s V 


where (f) are the Cartesian velocity components, v is the absolute fluid velocity, q the diffusive flux and S 0 are 
the volume sources. V is the computational cell volume, S are bounding cell surfaces, and v g is the grid velocity. If 
the grid is moving with time for the fluid-structure interaction problem, a space conservation law (SCL) is enforced 
during the grid deformation. 


— j dV = j v e ds 
dt v g 


(9) 


Related to the present liquid free surface problem, CFD-ACE+ contains a Volume of Fluid (VOF) module which 
is designed for applications involving two immiscible fluids. In the current application, the first fluid is LH2, LOX, 
or water, and the second fluid is gaseous H2, 02, He or water vapor, respectively. In the VOF module, a single set 
of momentum and continuity equations is solved, but different property sets are defined for each fluid. The volume 
fraction of one phase (in this case the liquid phase) is tracked throughout the solution to determine which fluid 
occupies each computational cell at any given time. In cells containing both fluids, a special routine is used to 
locate the shape, location and normal of the interface. When surface tension force is significant, its effect is applied 
in a conservative form. For time dependent simulations such as the present tank under oscillatory force, a special 
second order algorithm is used to update the volume fraction in a cell from one time step to the next. A second 
order geometric reconstruction scheme for the interface representation is employed to track the interface. This 
unique reconstruction algorithm is currently available only for structured quadrilateral and hexahedral grids. It is 
due to this algorithm that the present CFD solution is capable of capturing and maintaining a sharp interface 
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between the phases for very long period of time. The time step size is determined by the local Courant-Friedrichs- 
Lewy (CFL) number. Some of the validation studies and applications to space vehicle vertical and side sloshing can 
be found in references [9-11]. 

III. Validation Study for Smooth Wall Slosh Damping 

We will use the experimental data correlated into the industry standard as our baseline validation. Figure 1 
shows the experimental points we selected and the relationship to the correlation for the straight cylinders. The 
cylinder rdius varies from 3.8cm to 45.7cm and the contained fluid is water under normal gravity (lg) condition. 
Due to the variation of the cylinder size, the damping ratio from experiments varies from 0.334% to 0.0776%. Our 
goal is to determine the capability of the CFD solution in reproducing the damping trend and the accuracy of the 
prediction. 



Data 

point 

Tank radius, 
cm 

Kinematic 

viscosity, 

stokes 

Test liquid, 
[ reference] 

1 

15.2 

a 836 

SAE 10 -W oil [4.13] 

2 

3.8 

a 00929 

Water [4.4] 

i 

15.2 

U.U72T 

iserosene [«.«] 
Water [4.4] 

4 

7.6 

0.00929 

? 

15.2 

0. 00929 

Water 14. 13] 

6 

38.1 

0.0223 

Kerosene |4.13| 

7 

38.1 

0.00929 

Water [4. 13] 
Water [4.131 

8 

45.7 

0.00929 

4 

154.0 

0.UIN2V 

Unpublished lonvairr 
Astronautics data 


Figure 1. Experimental data and empirical correlation used for the current validation study 

Simulation Model 

The simulation model along with the computational grid are shown in Figure 2. Flere the first data point (#2) of 
a cylindrical tank with a radius of R=1.5” (3.8cm) and height of H=12” is modeled. The tank has no internal 
obstruction and all the walls are smooth. The top of the tank is open to ambient and is set as fixed pressure 
boundary condition. All other boundaries are set as non-slip walls. Shown in Figure 2 is the coarsest mesh with a 
total of 40,000 hexahedral cells. There are 60 cells covering the whole tank diameter. This ensures that there will 
be at least 120 cells to resolve the first natural mode and 60 cells to resolve the second mode. A symmetry condition 
is assumed so that only half of the cylinder is considered. Butterfly grid methodology is used to avoid the singularity 
at tank center line, to give better uniformity in grid spacing, and to ensure resolution of wall shear stress. One can 
see from Figure 2 the grid is packed near the wall to resolve the boundary layer which contributes to slosh damping. 

The initial liquid fill level is 2R, therefore the simplified expression (4) can be used. To determine the slosh 
damping, an initial inclined surface is specified with a slope of 5.7 deg. The time step size is controlled by the CFL 
number, which is set at 0.1. The maximum allowed time step size is 0.5 ms. The temporal scheme is 2 nd order 
Crank-Nicolson, and the spatial scheme starts with 2 nd order central with 1% damping (meaning a hybrid with 99% 
of 2 nd order scheme and 1% the 1 st order upwind). As discussed in the following section, for finer mesh the solution 
instability forces one to use only the first order spatial scheme in velocity, although the solution of Volume-of-Fluid 
is still second order accurate in space. 

3.8cm Radius Tank with Water 

A typical time dependent mass center coordinate of the liquid mass inside the tank, computed by integrating the 
Volume of the Fluid value in each cell, is shown in Figure 3. As the tank is stationary and there is no energy input to 
the system, the liquid surface experiences free oscillation. The amplitude of successive slosh oscillations 
represented by the total mass displacement deceases because of energy dissipation (and due to numerical dissipation 
inherent in the CFD solution). The damping ratio can be computed as: 
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( 10 ) 


y = — In 
2 K 



Where A is the mass center amplitude at i th cycle. The extracted slosh damping is 0.0109 or 1.09%, which is 3 
times higher than the experimental data of 0.00334 (0.334%) as seen from Figure 3. Indeed, it is very challenging to 
accurately predict the slosh damping using a CFD tool. This result indicates the requirement for grid refinement, 
which should reduce the numerical damping. 




Figure 3. Total mass center with time for 3.8cm radius cylinder and the predicted damping in comparison to 
experimental data. 


Figure 4 shows a series of grids that are used for the grid refinement study. The total number of cells are 40K, 
256K, 1 Million and 4 Million. As a results, the capability in resolving the first natural mode has increased from 
120 cells, to 200 cells, to 320 cells, and to 520 cells. 


The mass center response to the initial inclined surface for different grids of 40K, 256K and 1M are shown in 
Figure 5. The slosh frequency is seen to be almost independent of grid resolution, as slosh frequency represents the 
bulk fluid motion. However, the finer mesh of 256K gives less damping compared to the case of 40K cells. The 
predicted damping values for different grids are listed in Table 1. One can clearly see that there is a significant 
accuracy increase in the predicted damping value in comparison to test data with using the finer mesh of 256K: 
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0.312% vs. 0.334% (experimental data) (see Figure 5). The predicted damping value is nearly grid independent for 
1M cells: 0.34% vs. 0.312% (256K). The CFD simulation captured the fundamental damping physics: viscous shear 
near wall and dissipation in fluid motion inside the tank. The present grid refinement study suggests that with an 
increase in tank size, a finer mesh will be required. 


Increasing spatial resolution 

> 



Figure 4. Models for computational grid refinement study. 




Figure 5. Effect of grid resolution on the slosh mass response and predicted slosh damping for a cylindrical 
tank of 3.8cm (1.5”) radius. 


Table 1 Grid Resolution Effects on Predicted Slosh Damping 
for 3.8 cm Radius Tank 



40K 

256K 

1M 

Test 

Damping 

1.09% 

0.312% 

0.340% 

0.334% 


15.2cm Radius Tank with Water 

Next we consider a tank that is 4x larger, i.e. R=15.2cm (6”). We use the same grid and model as before, except 
that the geometry is scaled by 4x. 

Given in Figure 6 is the mass center response to the initial inclined interface for different grids of 40K, 256K and 
1M for this larger size tank. Due to the increase in tank radius, the slosh frequency has dropped. Here one observes 
the same trend: higher numerical damping for the coarse mesh of 40K and nearly grid independent damping from 
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256K and 1M cells. The predicted damping values for different grids are listed in Table 2. One notices almost 4 
times higher damping from 40K cells simulation indicating high numerical damping from the solution. The 
predictions using 256K and 1M cells are very close to each other and to the experimental value. Based on the two 
sizes of tank, the correlation of equation (4) has been reproduced. This validation study builds confidence that CFD 
simulation can indeed well predict the small values of slosh damping. 



Figure 6. Effect of grid resolution on the slosh mass response and predicted slosh damping for a cylindrical 
tank of 15.2cm (6”) radius. 


Table 2 Grid Resolution Effects on Predicted Slosh 
Damping for 14.5 cm Radius Tank 



40 K 

256K 

1M 

Test 

Damping 

0.895% 

0.206% 

0.210% 

0.199% 


45.7cm Radius Tank with Water 

The next tank has a radius of 45.7cm or 18”. It is 12x larger than the first tank (3.8cm). During our study we 
found that as tank size increases, the second order spatial scheme has difficulty in convergence due to extremely low 
physical dissipation in the system, and one is forced to use the first order upwind. It is should be pointed out that 
even though the spatial scheme for velocity is 1 st order, the scheme for Volume -of-Fluid which tracks the gas-liquid 
interface is still second order accurate. As a result of the use of 1 st order upwind, the predicted damping values 
from 256K or 1M cells are higher than the experiment data as observed from Figure 7. However, with the use of 
increasing cell numbers from 40K, to 256K, to 1M, and to 4M, one can clearly see the reduction in the predicted 
damping. With 4M cells the predicted damping value is very close to the experimental data. The call for increased 
mesh density shows the need for a highly scalable CFD code which is under development. Table 3 lists the 
predicted damping value and comparison with the experimental data. 
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Figure 7. Effect of grid resolution on the slosh mass response and predicted slosh damping for a cylindrical 
tank of 45.7cm radius 


Table 3 Grid Resolution Effects on Predicted Slosh Damping for 

45.7cm Radius Tank 



40K 

256K 

1M 

4M 

Test 

Damping 

0.438% 

0.162% 

0.140% 

0.0817% 

0.0776% 


IV. Validation Study for Slosh Damping with Ring Baffle 

Based on the above validation study, one can see that damping caused solely by viscous effects in a smooth tank 
is quite small; especially for even a moderately large tank. Values of y are generally equal to 0.001 or less when the 
tank diameter is 1 meter or so. Guidance, Navigation and Control analysis usually requires y to be considerably 
larger than 0.01 for stable flight. Consequently, some method of increasing the damping is required. For 
axisymmetric tanks, a common method is to attach a series of ring baffles to the tank walls (see Figure 8). The 
baffles provide a substantial degree of damping when the free surface is near one of the baffles. The spacing 
between the baffles is chosen so that the damping exceeds the requirements, regardless of fill level. 

Fuel-slosh damping by ring baffles in cylindrical tanks has been investigated extensively in the 1960’s [13-18], 
both theoretically and experimentally. A survey of damping measurements obtained in various experiments shows 
many apparent discrepancies. With the above validation, the purpose of this study is to present a CFD approach for 
fuel-slosh damping which brings all of the results together and provides a logical basis for prediction of damping in 
full-scale tanks. 
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Figure 8. Schematic of ring baffle and damping in a cylindrical tank 

The most widely used damping equation at present is the one obtained by Miles [17] which is based on 
experiments of Keulegan and Carpenter [18]. The equation is written as : 


83e -4 ' 60d/R a (3/2) 


' 


ct — Ag/ Ag 


2 (w/R) - ( w/R) ‘ 


(ID 


Here, d, w, A B , and a respectively denote the baffle depth, width, area, blockage ratio, while q, R, and AT 
denote the slosh wave amplitude, local tank radius, and tank cross-sectional area. O'Neill suggested that the wall 
side force amplitude (F) is easier to measure experimentally than q , and is linearly related the slosh wave amplitude 
in dimensionless form. Thus the following Miles-O'Neill equation is a convenient modification [14] . 

F = F/(pgR3) * 1.71 (q/R) -> £ = 2.16W 4 ' 60d/R a (3/2) ^ (12) 


Here, F , p, and g are the dimensionless side force, the liquid density, and acceleration of gravity. Since F is 
also easier to measure than q in CFD simulations as well, the above Miles-O'Neill form will be used to compare 
CFD results. 

Computational Model 

Two experimental sets were selected for comparison: the investigation of Silveira, Stephens, and Leonard [13] 
due to their measurement of slosh frequency shift, and that of O'Neill [14] due to his measurement of damping vs. 
dimensionless side force. Both sets include studies of slosh wave decay in cylindrical tanks and are outlined in [15]. 

A grid with a single baffle located at dimensionless height (h-d)/R = 1.25 was initially used to develop the 
simulation. Subsequently, two grids were created with a single baffle located at (h-d)/R = 2.0, where d is the 
distance between baffle and the quiescent free surface, and h is the free surface height (see Figure 10). This 
increased baffle height limits effects due to the tank bottom and is more consistent with experiments. The first grid 
has a baffle width ratio, w/R = 0.125 (a = 0.234, see equation (11)) while the second grid has w/R = 0.240 (a = 
0.422). 

Our previous grid refinement study found that at least four cells are needed along both the thickness and width of 
the baffle. To this end, a baffle thickness ratio of t/R = 1.5% was employed. While this is more than an order-of- 
magnitude thicker than that usually used in practice (t/R ~ 0.1%), the baffle thickness is generally not considered a 
critical parameter [16] and {t/R, t/w] « 1. The hyperbolic tangent edge feature of the CFD-GEOM program was 
used to blend the grid spacing defined by the baffle thickness up to that defined by the top and bottom faces of the 
cylinder. These considerations led to grids of approximately 60K cells. 

Since ideal transverse slosh occurs only along a single axis (x), half-cylindrical domains were selected in order to 
reduce the required computation time. The boundary conditions employed are no-slip along the tank walls and baffle 
edges, and symmetry along the xz plane. The w/R = 0.125 grid is illustrated in Figure 9, while Figures 10 and 11 
illustrate fluid motion in a typical simulation. The initial free surface with zero velocity is shown in Figure 10, 
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followed by its progression along with velocity vectors at increments of approximately 0.1 sec (along with velocity 
vector). In Figure 10, the liquid appears to be near a state of minimum velocity and maximum potential energy. 

□ 

Figure 9. Computational Fluid Dynamics Model for a Tank with Internal Ring Baffle 




Figure 10. Progression of Free Surface and Fluid Velocity (R=0.5ft, w/R = 0.125, d/R = 0.25, t = 0.0-0.3sec). 
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Figure 11. Progression of Free Surface and Fluid Velocity (R=0.5ft, w/R = 0.125, d/R = 0.25, t = 0.4-0.7sec). 

Baseline Results 


Simulations were first run with R=1.0m, (h-d)/R = 1.25, w/R = 0.15, and h/R = {1.500, 1.625, 1.750, 1.875, 
2.000} corresponding to d/R = {0.250, 0.375, 0.500, 0.625, 0.750}. The initial air-water free surface was specified 
as before. Values of damping ratio were obtained as follows by fitting an exponential curve to five positive peaks of 
the fluid center-of-gravity position, x cg and using the decay constant, p. 

x cg ,max = Ae- pt = Ae- ?Bt -> ^ = p/co (13) 

Center-of-gravity position data for d/R = 0.50 is shown in Figure 12. Results of damping ratio for all depth ratios 
are also shown in Figure 12 along with those predicted by the Miles-O'Neill Equation using the dimensionless side 
force of the middle (third) peak. As one can see the present CFD tool predicts the damping trend and value with 
baffle depth ratio fairly well. 



Figure 12. Computed force decay and damping for tank with baffle, and comparison with Mills-O’Neil 


In the next simulation set, the number of cells was uniformly increased from 50K to 75K and simulations were 
run for d/R = {0.250, 0.500, 0.750} only. The refined grid resulted in slightly lower values of damping ratio at each 
baffle depth ratio (2.7% max reduction). Unlike the smooth tank, slosh damping due to the baffle is less sensitive to 
the grid resolution. 
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For most simulations, the initial free surface slope was set to 0.2. Whenever possible, the simulation was allowed 
to progress for 2.5 cycles before using eg position data so that the liquid could settle into a more normal sloshing 
motion than afforded by the linear-surface initial condition. To observe the effect of tank radius, simulations were 

first run with the w/R = 0.125 grid scaled to R = {0.25, 0.50, 1.00} ft with d/R = 0.25. Since F can change rapidly 
at low depth ratios, values of the damping ratio are best compared with the Miles-O'Neill Equation by using the 
logarithmic decrement, 5, with n=2 as follows: 


5 = 2ni; = (l/n) In ( Ft / F i+n ) -> & = { 1/ ( 411) } In ( F i-\ / F i+ i ) (13) 

where n is the number of cycle observed. Thus, the Miles-O'Neill damping ratio at the i th peak of x cg can be 
compared to the CFD value calculated from the i-1 and i+1 peaks. Results are shown in Figure 13 for R = 

{0.25, 0.50 }ft. Results for R=1.00 ft are visually identical to those for R=0.5 ft indicating that baffle damping is 
independent of Re for tanks of practical size, as suggested by the equation. Again one can see fairly good agreement 
with Miles' equation from the present CFD simulations. 


R = 0.25 ft, w/R = 0.125, d/R = 0.25 



R = 0.5 ft, w/R = 0.125, d/R = 0.25 



(a) TankR=0.25ft (b) TankR=0.5ft 

Figure 13. Comparison of damping value predicted by CFD and that from empirical correlation of Miles- 
O’Neil at tank with sizes of R=0.25ft and 0.5ft 


Effect of Baffle Location 

For the next set of tests, the w/R = 0.125 baffle grid was again used to simulate sloshing, scaled to R=0.5ft with 
six different baffle depths below the free surface: d/R = {0.166, 0.250, 0.375, 0.500, 0.625, 0.750}. Note that the 

smallest depth ratio is very close to that of the experimental data by O’Neil [14] (d/R = 0.169). Data at F =0.15 are 
conveniently calculated without extrapolation and shown in Figures 14 and 15 along with the Miles-O'Neill 
damping. These plots are presented in terms of 8 for easy comparison with [13]. Selected experimental values from 
the reference have been extracted from [2] using the software, UN-SCAN-IT. The values have been scaled by a 
factor of 1.24 and included in Figures 14 and 15. This factor was selected so that the experimental 8 for d/R = 0.285, 
w/R = 0.241 matches the intersection of the CFD and Miles-O'Neill curves at the same depth ratio for w/R = 0.240. 

Data from [13] must be scaled because no value of F or q is specified in the reference. It can be seen that when the 
baffle is close to the free surface it provides high damping. The present CFD simulation correctly predicts the trend 
and damping value in comparison to the experimental data. 
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0.35 



Figure 14. Effect of baffle depth ratio on the slosh damping for baffle width ratio w/R=0.125 



Figure 15. Effect of baffle depth ratio on the slosh damping for baffle width ratio w/R=0.24 


IV. Conclusion 

The present study represented a continuing effort in validating CFD technology in modeling space vehicle 
propellant sloshing dynamics. This effort has clearly illustrated the soundness of the CFD approach in modeling the 
detailed fluid dynamics of tank sloshing and the excellent accuracy in extracting slosh mechanical properties, 
especially sloshing damping for smooth walls for different tank configurations and baffled tanks. The present study 
shows that CFD technology is especially applicable to the future design of propellant tanks, as there is no previous 
experimental data available for the same size and configuration. The previous analytical solutions are generally 
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obtained for the linear analysis, where the amplitude of the slosh wave is small. CFD solution can go beyond the 
linear stability regime, and predict what will happen after surface breakup. CFD simulation also has no limitation 
on the gravitational acceleration, and can be applied to a realistic flight model. 

Application of the CFD modeling tool has great potential in the design of propellant tank with multiple baffles. 
In general, the practice to estimate the total damping for multiple baffles is by a linear superposition of the estimated 
damping contribution from each individual baffle. This approach is considered adequate so long as the baffle 
separation distances exceeds their widths (S > W). However, experimental data involving multiple baffles are 
limited. CFD simulation could be used to determine the validity and applicability of such rules. Finally, CFD 
simulation could also be used to help develop a relation between slosh frequency and baffle parameters. Using 
baffles to shift frequency could be an alternative to tank compartmentalization. 
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